Star-unitary transformations. 
Prom dynamics to irreversibility and stochastic behavior 
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We consider a simple model of a classical harmonic oscillator coupled to a field. In standard ap- 
proaches Langevin-type equations for bare particles are derived from Hamiltonian dynamics. These 
equations contain memory terms and are time-reversal invariant. In contrast the phenomenolog- 
ical Langevin equations have no memory terms (they are Markovian equations) and give a time 
evolution split in two branches (semigroups), each of which breaks time symmetry. A standard 
approach to bridge dynamics with phenomenology is to consider the Markovian approximation of 
the former. In this paper we present a formulation in terms of dressed particles, which gives exact 
Markovian equations. We formulate dressed particles for Poincare nonintegrable systems, through 
an invertible transformation operator A introduced by Prigogine and collaborators. A is obtained 
by an extension of the canonical (unitary) transformation operator U that eliminates interactions 
for integrable systems. Our extension is based on the removal of divergences due to Poincare res- 
onances, which breaks time-symmetry. The unitarity of U is extended to "star-unitarity" for A. 
We show that A-transformed variables have the same time evolution as stochastic variables obey- 
ing Langevin equations, and that A-transformed distribution functions satisfy exact Fokker-Planck 
equations. The effects of Gaussian white noise are obtained by the non-distributive property of 
A with respect to products of dynamical variables. Therefore our method leads to a direct link 
between dynamics of Poincare nonintegrable systems, probability and stochasticity. 

PACS numbers: 02.50.Fz, 05.40.-a, 05.70.Ln 
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I. INTRODUCTION 

In classical physics the basic laws are time reversible. If 
we know the Hamiltonian, then we get Hamilton's equa- 
tions of motion which describe the time evolution of the 
system in a time reversible, deterministic way. On the 
other hand, we see time irreversibility and stochastic be- 
havior everywhere. How to bridge the gap between the- 
ory and reality has been the subject of many discussions. 

The main problem is how to extract irreversibility and 
stochasticity out of Hamilton's equations of motion. This 
will be the subject of this paper. Our approach is an 
extension of canonical transformations to define dressed 
particles or quasiparticles Q. 

We consider Hamiltonians that can be written as 



H = H + XV. 



(1) 



The first therm Hq describes a set of noninteracting 
"bare" units while the second XV describes their interac- 
tions (A is a dimensionless coupling constant). Specif- 
ically, we will consider the one-dimensional Fricdrichs 
model Q , describing a classical harmonic oscillator (bare 
particle) coupled to an infinite set of bare field modes 
(heat bath). This model is closely related to the Caldeira- 
Leggett model Q, which has been extensivelyused to 
study quantum Brownian motion 0, 0, S S 0- • 



In general, bare particles follow a complicated motion, 
due to their interactions. In order gain a physical insight 
into their behavior, and also to simplify the equations of 
motion, one can introduce a change of phase-space vari- 
ables (a canonical transformation). The new variables 
describe renormalized entities, or quasiparticles. After 
solving the equations for quasiparticles, one may apply 
the inverse canonical transformation to get the solutions 
of the original equations of motion. For the Friedrichs 
model the quasiparticle consists of the original particle 
surrounded by a "dressing" cloud of field modes. 

For integrable systems, one can construct transforma- 
tions U that completely eliminate the interactions. They 
bring us to a description in terms of free quasiparticles. 
This is readily seen starting with the Liouville equation 



.0 



(2) 



where Ljj = i{H, } is the Poisson bracket with the 
Hamiltonian. Similar to Eq. Q the Liouvillian is written 
as a free term plus interaction, Lh = Lq+XLv- Applying 
U on both sides of the Liouville equation we get 



i^-Up = ULnU-'Up 
at 



l dl P 



(3) 
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where 

p=Up, Lo^ULhU- 1 . (4) 

The transformation U is constructed in such a way that 
Lq has the same form as the non-interacting Liouvillian, 
with renormalized frequencies. Eq. Q gives the time 
evolution of the free dressed particles. 

If U can be constructed through a perturbation ex- 
pansion in A, we say the system is integrable in the sense 
of Poincare. For these systems we can keep a one-to- 
one correspondence between the original variables and 
the transformed variables. With a suitably defined inner 
product between dynamical variables and ensembles, we 
can define the hermitian conjugate transformation U'. 
One then finds that U is unitary: U' = U^ 1 . The trans- 
formation U thus preserves the time-reversibility of the 
original Liouville equation. 

Now, if all systems were integrable in Poincare's sense, 
this would mean that all the phenomena we observe in 
nature are equivalent to free motion. This would be hard 
to reconcile with the existence of dissipative phenomena, 
which are essential for the appearance of bifurcations and 
self-organization E3 • However for most systems one can- 
not construct U by perturbation expansions, due to the 
appearance of resonances. Resonances give vanishing de- 
nominators leading to divergences. These divergences 
were discovered by Poincare, so we will refer to them 
as Poincare divergences [hereafter, whenever we speak 
of integrability or nonintegrability, it will be meant in 
Poincare's sense]. 

It is precisely for Poincare's nonintegrable systems 
that we see irreversible and stochastic behavior, such 
as Brownian motion. One of the main developments 
of the Brussels-Austin groups led by I. Prigogine has 
been to show that one can systematically eliminate the 
Poincare divergences by regularization of denominators 
[H 13 El El III As a result of this regularization 
time-symmetry is broken and one obtains a new type of 
transformation A that replaces U. This gives a quasi- 
particle description leading to stochastic or kinetic equa- 
tions, such as the classical Langevin or Fokker-Planck 
equations, respectively. To see this, we operate A on the 
Liouville equation 

d 

i—Ap = KL H \- l kp 
at 

d 

= °P ( 5 ) 

where 

p = Ap, = ALhA -1 . (6) 

8 is now a collision operator as used in kinetic theory. 
The A transformation gives a probabilistic description, 
which is irreducible to trajectories in classical mechanics 
or wave functions in quantum mechanics. If we integrate 
out the field variables, Eq. (J5J becomes, e.g., an exact 
Fokker-Planck operator. Through the A transformation 



we can also describe dressed unstable states in quantum 
mechanics 0] . One can define as well an Tt- function 
that has strict monotonic behavior [Tlj . 

For the Friedrichs model we have both integrable and 
nonintegrable cases, depending on whether the spectrum 
of the field modes is discrete (finite volume L with pe- 
riodic boundaries) or continuous (L — > oo). In the first 
case we have cyclic (although complicated) motion of the 
particle, as the field comes back to the particle through 
the periodic boundaries. In the second case, the field does 
not come back. A Poincare resonance emerges, since the 
energy of the particle is embedded inside the continuous 
spectrum. The emission of the field from the particle 
leads to radiation damping. Conversely, the particle is 
excited when it absorbs the field. 

To understand the breaking of time symmetry, we note 
that when there are Poincare resonances, i.e. in the limit 
L — ► oo, the solutions of Hamilton's equations for the 
bare particle contain a dominant decaying (Markovian) 
component oriented either towards the future or the past, 
or both, depending on the initial conditions. Taking the 
well-known Markovian approximation, one finds that the 
equation of motion for the bare particle is split into two 
branches, one for t > and another for t < 0, corre- 
sponding to two semigroups. As a whole the time rever- 
sal invariance of the motion is kept, but if we pick cither 
branch, time symmetry is broken. 

To obtain this splitting into two semigroups for the 
bare particle we have to make approximations. In con- 
trast, in terms of the dressed particle defined through 
A, this is an exact property. The analytic continuation 
of U can be made to either the upper or lower complex 
frequency planes, giving exact Markovian equations that 
generate the t < or t > semigroups, respectively. 
Once we fix the analytic continuation, time symmetry is 
broken. 

The A-transformed functions involve generalized func- 
tions, or distributions (examples are the "Gamow modes" 
presented in Sec. 0. If the initial unperturbed functions 
formed a Hilbert space, the transformed functions are no 
more in this Hilbert space. In its transformed domain 
Lh behaves as the dissi pat ive collision operator 9 with 
complex eigenvalues fl3l Il4j| . 

In contrast to U , A is no more unitary. Instead, it is 
"star unitary" [Tl|,E3- Furthermore, while U is distribu- 
tive with respect to multilplication, A is non-distributive. 
As we will see, these properties allow us to describe 
damping and fluctuations associated with noise. 

A basic requirement on A is that it is invertible. This is 
connected with the star unitarity of this transformation 
[see comments below Eq. (|77ll ]. In addition to this, our 
construction of A is based on the following requirements: 

(1) The A transformation is obtained by analytic contin- 
uation of the unitary transformation U . When there are 
no resonances, A reduces to U . 

(2) A preserves the measure of the phase space. 

(3) A maps real variables to real variables. 

(4) A is analytic with respect to the coupling constant A 
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at A = 0. 

(5) A leads to closed Markovian kinetic equations. 

We will focus on the dynamical observables of the par- 
ticle. The action of A will be restricted to the subset of 
phase space functions depending only on the particle de- 
grees of freedom. Within this subset we obtain an exact 
and invertible A transformation. 

Furthermore, we will consider the thermodynamic 
limit of the field modes, where the total energy of the 
field is an extensive variable. Then the average action 
(Jfc) of each field mode k satisfies 0] 



(J k )-O(L ) 



(7) 



for L — > oo. The total energy of the field 

E f = u k (J k ) -> A f dkio k (J k ) ~ 0(L) (8) 
k 

is proportional to the volume L. [This does not neces- 
sarily imply that the field is Gibbsian]. The existence 
of the thermodynamic limit requires an initially random 
distribution of the phases of the field modes . 

A different situation occurs if the total energy of the 
field is in a non-extensive variable. Then we have (Jk) ~ 
i.e., we have a vanishing energy density. We will 
not consider this case in this paper. 

In the extensive case, in addition to the damped os- 
cillation, the particle undergoes an erratic motion due 
to the excitation caused by the field. This erratic mo- 
tion includes a Brownian motion component, which is 
Markovian. The initial randomness of the phases of the 
field modes is a necessary condition for the appearance 
of Brownian motion. In addition it is essential that the 
field resonates with the particle. We need Poincare res- 
onances. Under these conditions A permits us to isolate 
the damping and the Brownian component of the motion. 

Our approach can also be formulated in terms of com- 
plete sets of projection operators II'"), that permit to 
decompose dynamics into a set of orthogonal "subdy- 
namics" (see Sec. IV1|) . Essentially, we introduce a gen- 
eralized basis that permits us to analyze the motion in 
terms of strictly Markovian components. In our case, we 
study the component that describes Brownian motion. 
The other components, clumped together, give what is 
usually called non-Markovian (memory) effects fl4] . The 
Brownian component is independent of the initial corre- 
lations between the particle and the bath, and in this 
sense, it has a "universal" character. 

The results presented here are based on Refs. [15|, 
where we constructed A for the quantum Friedrichs 
model. The main subject in these papers was the de- 
cay of unstable particle states. We showed that the 
A transformation permits us to isolate the exponential 
(Markovian) component of the decay, which occurs when 
the energy of the field is non-extensive. The remaining 



(non-Markovian) component gives the Zeno effect 
and long tails 20], which are conneted to the appear- 
ance of a dressing cloud around the bare particle. The 



dressed unstable state defined through A has a real aver- 
age energy and gives an uncertainty relation between the 
lifetime and energy (see also Similar considerations 

can be applied in classical mechanics (2^| . 

The present paper is organized as follows. In Sees, [n] 
- IIHI we introduce the Friedrichs model and we discuss 
the equations of motion of the bare particle. In the sub- 
sequent Sections we study the evolution of renormalized 
(dressed) particle variables. We consider first (Sec. IIV|) 
the integrable case where the spectrum of the field is 
discrete. We can then define the renormalized variables 
through the unitary transformation U. In the continu- 
ous spectrum limit the system becomes nonintegrable. In 
Sec. as a first step to introduce A, we extend the renor- 
malized particle modes in the discrete case to the decay- 
ing "Gamow" modes in the nonintegrable case. In Sec. 
^Jwe construct A. In Sec. IVIII we show the correspon- 
dence between the solution of Langevin equation (with 
Gaussian white noise) and the A transformed variables. 
Finally, in Sec. IVIIII we derive a Fokker-Planck equation 
for the A transformed distribution function. Details of 
calculations are given in the Appendices. 



II. THE CLASSICAL FRIEDRICHS MODEL 

We consider a classical system consisting of a har- 
monic oscillator coupled to a classical scalar field in one- 
dimensional space. A quantum version of this model has 
been studied by Friedrichs 0, among others. 

We write the Hamiltonian of the system in terms of 
the bare oscillator and field modes ai and a^, 



\ " 



^uj k a* k a k + V k (a\a k + ajajl), 

k k 



(9) 



with a given constant frequency u)\ > for the harmonic 
oscillator (particle), c = 1 for the speed of light, and 
u>k = \k\ for the field. When A is small we can treat 
the interaction potential as a perturbation. We assume 
the system is in a one-dimensional box of size L with 
periodic boundary conditions. Then the spectrum of the 
field is discrete, i.e., k = 2irj /L where j is an integer. We 
assume that 



u>i =/= LOk , for all k 



(10) 



by 



The volume dependence of the interaction V k is given 



Vk 



(11) 



where Vk = 0(1). We assume that Vk is real and even: 
Vk — V-k- Furthermore, we assume that for small k 



1/2 



(12) 
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An example is the Drude-Ullersma form 0: 



1/2 



1 + 



(13) 



where lom is the cutoff frequency of the bath. 

To deal with the continuous spectrum of the field we 
take the limit L — > oo. In this limit we have 



2tt 



E 



dk, J±5 kfi ^5(k). (14) 
Ztt 



We will often use the summation sign with the under- 
standing that we replace it by an integral in the limit 
L — > oo. 

The bare modes a±, a k satisfy the Poisson bracket re- 
lation 



where 



i{a a , a%} = S af: 



iff ] = \^ \ d f dg d 9 df 
^ Yda r dat da r dat 



(15) 



(16) 



[the sum includes the discrete index r = 1 as well as the 
running index r = k] . The bare modes are related to the 
position x\ and the momentum pi of the particle as 



a i 



xi 



Pi 



mu)\ 



ipi 
mail ' 



(17) 



1 



= (ai + a[), 



(ax - a*) 



(18) 



and to the field (j>(x) and its conjugate field n(x) as 

1 xl/2 



E(^) ' (a k e ikx +ale- ikx ), (19) 



7r(x) 



/Wfcx 1/2 

"2L> 



(a k e 



ik.r 



(20) 



The field <j)(x) corresponds to the transverse vector poten- 
tial in electromagnetism, while 7r(x) corresponds to the 
transverse displacement field. Our Hamiltonian can be 
seen as a simplified version of a classical dipole molecule 
interacting with a classical radiation field in the dipole 
approximation |28j |. For simplicity we neglect the in- 
teractions proportional to a\a k and a^a^, which cor- 
respond to "virtual processes" in quantum mechanics. 
This approximation corresponds to the so-called rotat- 
ing wave approximation [J, |24| . If we incorporate the 
virtual processes, then we obtain the classical version of 
the Caldeira-Leggett model. 

We note that we have an cu k — uj- k degeneracy in 
our Hamiltonian. To avoid some complexity due to this 



degeneracy, we rewrite our Hamiltonian in terms of new 
bare modes as [2^ 

H = Wi^gi + ^2 u kq* k qk + A V k (q{q k + q x ql), (21) 

k k 

where 

q\=a\, q k 



_ / (a k + a-k)/V%, for k > 0, 
(o fc - a_ fc )/V2, forfc<0, 



(22) 



V, 



V2V k , 
0, 



V k . 



for k > 
for k < 0, 



(23) 



(24) 



In this form the mode q k with negative k argument is 
completely decoupled from the other degrees of freedom. 
The new bare modes also satisfy equations (fHfl) . JTBJl. 

In the subsequent sections we will use the following 
notations. We define action and angle variables J s ,ct s 
through the relation 

(25) 



q s = V^e M % s = 1, k 
We define T as the set of all modes, 
r = (qi,ql, ...,qk,q k , ■■■) 



(26) 



and T s = (q s ,qt), with s — l,k, as the set of particle or 
field modes. We will also denote as the set of all field 
modes Tf = {Tfc}. We use the notation dT for the phase 
space volume element and dTi, dTf for the particle and 
field components of dT, respectively 

dT = dTxdTf 
dTi = dJ\da\, 

We define as well 
5(T-V) 

= 8(Jx - J[)8(ai - a\ 



dV f = Y\_dJ k da k . 



(27) 



<Y[6(J k -J' k )5(a k -a k ). 



(28) 



We consider ensemble averages as inner products: 

(F) = ((F\p)) = J dTF(Typ(T). (29) 
For an operator O the Hermitian conjugate is defined by 



((F\Op)) = ((p\O^F)Y 



(30) 



As mentioned in the Introduction, for our model we can 
have both integrable (L finite) and nonintcgrable cases 
(L — > oo). In the first case there are no resonances (see 
Eq. (|10|) ~) and as we will see, the the system is integrable 
in the sense of Poincare. In the second case the system 
can become nonintegrable in Poincare's sense, due to the 
emergence of the resonance u)\ — w k between the fre- 
quencies of the particle and the field. This distinction is 
essential in our construction of dressed particle modes. 
Before coming to this, we will briefly consider the equa- 
tions of motion for the bare modes. 
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III. EQUATIONS OF MOTION OF THE BARE 
PARTICLE MODES 



The dynamical equations of motion of an oscilla- 
tor coupled to a field have been studied by many au- 
thors, mainly using the quantum Caldeira-Leggett model 
0, 0, EL 0- • Here we will write the equations for the 
Friedrichs model. In contrast to the phenomenological 
equations describing Brownian motion |25j . these equa- 
tions have memory terms (i.e. they have time-dependent 
damping and diffusion coefficients), and the time evolu- 
tion they generate forms a group, since they are equiva- 
lent to Hamiltonian dynamics. 



A. Non-Markovian Langevin equation 

Starting from the Hamiltonian equations 

q,(t) = exp(iL H t)q a {0) => q s (t) = iL H q s (t), (31) 
we can obtain the exact time evolution of the modes q r 



q s {t) = X)/«r(*)3r(°) 



(32) 



where f rs {t) are complex functions (see Appendix Ia"|) . 
We will focus our attention on the particle modes: 

qi{t) = fll(t)qi(0) + Y,.hk(t)q k (0), 

k 

qi{t) = fn(t)qi(0) + ^2fik(t)q k (0). (33) 

k 

Solving for qi(0) in the first equation and replacing the 
result in the second equation we get 



qi{t) = -izx{t)qx{t) + R(t) 



where 



R(t) 



«^ln fn{t), 
Y,h k (t)q k {0), 



h k (t) = fik(t)+izi{t)f lk {i). 



(34) 

(35) 
(36) 

(37) 



Eq. (|34|) is a non-Markovian equation, because of the 
time dependence of the coefficients. The function z\(t) = 
uj\{t) — ij(t) gives the instantaneous frequency Qi(t) and 
damping rate ^(t) of the oscillator [we note that damping 
appears only in the non-integrable case]. R(t) is an er- 
ratic function, since it depends on the initial states of all 
the field modes <?fc(0) (see Appendix [BJ). It plays the role 
of noise. In general, this is colored noise, as the function 
R(t) has memory in the auto-correlation, 



where 



(R*(t)R(t')) ^ fori 
means ensemble average. 



(38) 



B. Non-Markovian Fokker-Planck equation 

We can also derive a non-Markovian equation for the 
particle distribution function 



pi(Ti,t)= J dV }P {T,t) 



(39) 



which allows us to calculate averages of functions G(T{) 
depending only on the particle modes. We assume that 
G(Fi) is a smooth real function of Ti that vanishes at 
|<7i| = oo and is expandable in the infinite series 



G(q u qt)= J2J2G 



*m n 
mull Hi ■ 



(40) 



We assume as well that p(T, t) is factorized at t = into 
independent particle and field mode functions and that 
the field distributions depend only on the actions. In 
other words, we have 



p(r ) o) = Pl (r 1) o)IIftW. 



(41) 



Our final assumption is that the volume of the system is 
large, so we neglect terms of order 1/L. This approxi- 
mation becomes exact in the continuous spectrum limit 
L — > oo, i.e., in the nonintegrable case. We consider the 
extensive case discussed in the Introduction. 

Based on Eq. Q34[l we then obtain the non-Markovian 
equation (see Appendix Q 



iz^^-qx-izKt^qt (42) 
dqi dq{ 



where 



m = E 



—izl{t) + iz\{t) + 



d_ 

di 



\hk(t)\ 2 (J k ) . (43) 



Eq. 142|) is of the Fokker-Planck type, but with time de- 
pendent coefficients. 

The equations derived in this section are reminiscent 
of the phenomenological equations for Brownian motion. 
However, the phenomenological equations have quite im- 
portant differences: they are Markovian, they break 
time-symmetry and they describe stochastic processes. 
One can derive the phenomenological equations using 
approximations, such as the Markovian approximation. 
This is shown in Appendix iDl 

In the rest of the paper we will study the dynami- 
cal evolution of renormalized modes. In contrast to the 
bare modes, the renormalized modes obey exact equa- 
tions having the same evolution as the phenomenological 
equations. As a preparation, we first consider the inte- 
grable case. 
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IV. UNITARY TRANSFORMATION FOR 
INTEGRABLE CASE 

In this Section, we present the properties of the canon- 
ical transformation U that diagonalizes the Hamiltonian 
in the discrete spectrum case, when the size of the box 
L is finite. Later we will extend U to A through ana- 
lytic continuation, for L — > oo. In the integrable case we 
can find renormalized modes Q s , Q* that diagonalize the 
Hamiltonian through U . The new modes are related to 
the bare modes as 



Q s = IP q s for s = 1, k. 



(44) 



in one-to-one correspondence. The operator U is unitary 

u- l = W. 

The Hamiltonian is diagonalized as 



(45) 



where uj a are renormalized frequencies. 

The new modes satisfy the Poisson bracket relation 



i{Qr, Q* s } = <5r 



(46) 



Since the interaction is bilinear in the bare modes, the 
new modes can be found explicitly through a linear su- 
perposition of the bare modes j^. For the particle we 
obtain, from the equation i{H, Qi} = —Q\Q\i 



Qi=Nl /2 ( qi +\J2c k qk) 



where 



V k 



UJl - UJ k 



(47) 

(48) 
(49) 



The renormalized frequency uj\ is given by the root of 
the equation 



T)(ui) = 0, Tf(z) 



Z — UJ\ 



E 



z - to k > 



(50) 



that reduces to uj\ when A = 0. For the field modes one 
can also find explicit forms (see Appendix [SJ. 
The perturbation expansion of Eq. (|47|l yields 



h = uUi = qi + 



UJl — UJk 



-q k + 0(X 2 ) (51) 



When the spectrum is discrete, the denominator never 
vanishes; each term in the perturbation series is finite. 
This implies integrability in the sense of Poincare: U can 
be constructed by a perturbation series in powers of A™ 
with n > integer. In other words, U is analytic at 
A = 0. 



Since the transformation U is canonical, it is distribu- 
tive with respect to multiplication 

UU r q* s = [UUr}[tfq* s ]=QrQ* s (52) 
Hence we have 

uh = u[J2 io s Q:Q s ] = J2 Q ^q s = (53) 

The transformed Hamiltonian UH has the same form 
of the unperturbed Hamiltonian Hq, with renormalized 
frequencies. 

The canonical transformation can also be introduced 
on the level of statistical ensembles p, as shown in the 
Introduction. In Eq. J2J we have 



L p = iU{H, p} = i{UH, Up} 

Up 



dq s 



(54) 



where in the second equality we used Eq. I|53|l and the 
property of preservation of the Poisson bracket by canon- 
ical transformations [2(j . Hence the transformed Liouvil- 
lian L does not contain any interaction terms. Ensem- 
ble averages over this transformed density function p can 
thus be easily calculated. For example for 



i^- t ((x 1 \U\p))=i^(x 1 ) 



(55) 



and similarly for (pi) we get, after substituting Eq. I |18(l 
and integrating by parts, 



4 = (Pi) . 4 (pi) 
at m at 



-muj 1 [Xi) 



(56) 



These are the equations for the free harmonic oscillator 
(with renormalized frequency uj\ and renormalized mass 
rn = rawi/iDi). The interaction with the field is elimi- 
nated. 

Note that the normal modes are eigenfunctions of the 
Liouvillian Lq, 



L <7i = -uiqi, L q* = uiiq^. 
This leads to 

L H Q 1 = -uj 1 Qi, L H Q* 1 =ui 1 Q*. 
For products of modes we have 



L Q qt m qi 



[{m-n)u>x]ql m q^ 
[(m-n)wi]Oi- 



(57) 



(58) 



(59) 



Finally, we note that from distributive property Eq. 
(15211 we have 



uUl m qi = (f/ t g 1 * m )(t/ t gi)- 



(60) 
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V. NONINTEGRABLE CASE: GAMOW MODES 

Now we consider the continuous spectrum case, where 
the particle frequency w\ is inside the range of the con- 
tinuous spectrum u>k- In this case, by analytic contin- 
uation of Qi and Q\ we can get new modes which are 
eigenf unctions of the Liouvillian with complex eigenval- 
ues. These modes are called Gamow modes. Gamow 
states have been previously introduced in quantum me- 
chanics to study unstable states [2]j-[32. In classical 
mechanics, Gamow modes have been introduced in Ref. 
|22T |. In this Section we present the main properties of 
Gamow modes, which will be used for the construction 
of A. 

When we go to the continuous limit we restrict the 
strength of the coupling constant A so that 



ak < u>i, 



(61) 



Then the harmonic oscillator becomes unstable. In this 
case we have radiation damping. If Eq. (|fi If) is not sat- 
isfied, then we go outside the range of applicability of 
the "rotating wave approximation" (see comment after 
Eq. l]2Up'l as the Hamiltonian becomes not bounded from 
below, and gives no radiation damping [3^ |. 

In the continuous spectrum case, divergences appear 
in the construction of U, due to resonances. For exam- 
ple, the denominator in Eq. I|51|l may now vanish at the 
Poincare resonance u>i = ujk- We have a divergence in 
the perturbation expansion in A. To deal with this di- 
vergence, we regularize the denominator by adding an 
infinitesimal ±ie. Then we get 



Qi = gi + 5Z 



u>i — u>k ± ie 



q k + 0(\ 2 ). 



(62) 



In the continuous limit the summation goes to an inte- 
gral. We take the limit L — > oo first and e — > oo later. 
Then the denominator can be interpreted as a distribu- 
tion under the integration over k 



1 



1 



=F inS(wi - LOk) (63) 



where V means principal part. 

The introduction of ie in the continuous limit is re- 
lated to a change of the physical situation. In the dis- 
crete case the boundaries of the system cause periodicity 
in the motion of the particle and the field. In contrast, in 
the continuous case the boundaries play no role. In the 
continuous limit we can have damping of the particle, as 
the field emitted from the particle goes away and never 
comes back. And we can have Brownian motion, due 
to the interaction with the continuous set of field modes. 
The continuous limit may be well approximated by a dis- 
crete system during time scales much shorter than the 
time scale for which the field goes across the boundaries. 

In the continuous limit we can have damping of the 
particle either toward the future or toward the past. This 



corresponds to the existence of the two branches ±ie in 
Eq. (|62|) . Breaking of time symmetry is connected to 
resonances |3~H . 

As shown in |30j ]. continuing the perturbation expan- 
sion (|62|l to all orders one obtains new renormalizcd 
modes (Gamow modes) associated with the complex fre- 
quency 



Z\ = Wi — 27 



(64) 



or its complex conjugate z*. Here Q\ is the renormalized 
frequency of the particle, and 27 > is the damping rate. 
The complex frequencies are solutions of the equation 



. / ,//. x2v l =0 . 



Tj {u>) = u) — u)± — / dk 



(z - w k )l 



(65) 



The + (— ) superscript indicates that the propagator is 
first evaluated on the upper (lower) half plane of z and 
then analytically continued to z = u). 

The new modes for the — ie branch in Eq. I|62|l are given 

by 



Qi = Nl /2 [ qi + Xj^CkQk], 



(66) 



Ck = 



V k 



AT 1 = (1 + A 2 ^4)- 1 (67) 



(Z - bJk)ii 

and its complex conjugate, satisfying 

L H Qx = -zxQx, L B Ql = zlQ\. (68) 
The mode Q\ decays for t > as 



e iLntQ* _ e i*itQ* = g(«5i-7)*Q* 

(and similarly Q\). 



(69) 



The modes for the +ie branch are given by 

QI=^ /a [gJ+A£>9*] (70) 



and its complex conjugate, satisfying 

LhQ* = Z1Q1, LhQi = —z*Qi. 



(71) 



These modes decay for t < 0. 

The modes we have introduced have quite different 
properties from the usual canonical variables. Their Pois- 
son brackets vanish 



i{QuQl} = i{Qi,Qt} = 0. 



(72) 



However the modes Q\ and Q\ are duals; they form a 
generalized canonical pair 



i{Qi,Ql} = l. 



(73) 



This algebra corresponds to an extension of the usual Lie 
algebra including dissipation. An analogue of this algebra 
has been previously studied in quantum mechanics p?j- 
321, in terms of non-Hilbertian bras and kets. 
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VI. THE A TRANSFORMATION 

Using the above results we now introduce A. In this 
paper we will restrict the action of A to products of par- 
ticle modes of the form q* m g™. This will be enough to 
calculate renormalized functions of the particle variables 
(expandable in monomials), which will lead us to the 
Langevin and Fokker-Planck equations. The action of A 
on more general functions, including field modes will be 
considered elsewhere (see also [T5I l22| ~) . 



A. Defining A through its action on particle modes 

First recall that in the integrable case the renormalized 
particle modes are related to the original modes as 



qi, 



Ql = W 



Qi 



(74) 



For products of modes we have as well the relation (|6U|I . 
In the continuous spectrum limit we come to the nonin- 
tegrable case. As seen in the previous Section, we elimi- 
nate Poincare divergences in single renormalized particle 
modes by analytic continuation of frequencies to the com- 
plex plane (i.e. u>\ goes to z±) leading to Gamow modes. 
There are two branches for the continuation, namely 



Qi 
Qi 



(75) 



Corresponding to these extensions, we introduce A, the 
extension of U in Eq. i|74|) . 



Qi 

Qi 



A" 1 



Qi, 



= A- a 



Qi- 



(76) 



These relations partially define A, by its action on single 
particle modes (a more complete definition is given be- 
low). This definition satisfies the requirements (1), (3) 
and (4) given in the Introduction. We will comment on 
the remaining requirements (2) and (5) below. Note that 
A' ^ A -1 is not unitary. Instead, it is "star-unitary," 



A" 1 = A*. 



(77) 



In our case, where we restrict the action of A to particle 
modes, star conjugation has a simple meaning. It simply 
means taking hermitian conjugation and changing ie => 
— ie, so we have, e.g., [A*(ie)]qi = [A^(— ie)W For the 
general definition of star conjugation, see |3(j and 
references therein. 

Due to star-unitarity, the existence of the star- 
conjugate transformation A* guarantees the existence of 
the inverse A -1 . 

As mentioned above we are interested not only in the 
renormalized modes, but also the renormalized products 
of modes, 



At 



?r?i. a 



-YrQi- 



(78) 



For the integrable case, renormalized products of 
modes can be easily calculated thanks to the distribu- 
tive property (|6U[) . However, as shown below, for the 
nonintegrable case products of Gamow modes give new 
Poincare divergences. Hence, due to the requirement (4) 
stated in the Introduction, the A transformation has to 
be non-distributive. This means that we still have to 
define the action of A on products of particle modes. 

Let us first consider the transformed product A^q*qi. 
Later we will generalize this to obtain the expressions 
(|7S|) . If A* were distributive, A^q*qi could be expressed 
as the product Q\Q \ — (A^ ql)(A^ qx) . However, as we 
show now, this expression gives Poincare divergences in 
the thermodynamic limit. We have 

QlQi 

= {NiM + X^clqDiqx + X^CkQk) 



\Ni\(qlqi + Xql c k q k + Xqi c 



lit 



+ A 2 ^4 Cfc ,^+A 2 ^| Cfc | 2 ^ fc )- (79) 

where the prime in the summation means k ^ k' . Going 
to the continuous limit and taking the ensemble average 
with an ensemble p the last term becomes 



E 



Cfc| {QkQk) 



dk 



(z - ujk)ii 



(J k ) (80) 



where (Jk) = ((qlq k \p)). This term has a non-vanishing 
finite value in the limit L — ► 00 if Eq. Q is satisfied. 
Furthermore, if the ensemble p belongs to the class of 
ensembles with ^-function singularities in the wave num- 
ber k, then Eq. ||SD) is non-negligible as compared to the 
average of the q^Qk' term in Eq. (|79H . For this class of 
ensembles the point contribution k = k' is as important 
as the integration over k' [la. 1 1 ll l35j : 



Eto')-U)-0(L°). 



(81) 



(see Appendix[EJ . This type of ensembles with 6- function 
singularities is by no means atypical. An example of 
this class of ensembles is the Gibbs distribution. For 
ensembles in this class, we have well defined intensive 
and extensive variables in the thermodynamic limit |35j . 
To lowest order we have in Eq. I|8U|I . 



Xv k 



Xv k 



(z - Wfc)J, 
which leads to 
Xv k 2 



uji — u k + ie 



(z - 



\uji -u k + ie\ 



0(X 3 



0(A 4 ) 



(82) 



(83) 



X v k 8(uji - uj k ) + 0(A ) 



CO. 
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This diverges when e — > 0. Hence Eq. I|80[) is nonana- 
lytic at A = due to the resonance at uji — u> k . We 
have Poincare divergence in the perturbation series of 
(Attf)(At ?1 ). 

We note that when the energy of the field is non- 
extensive, we have (J k ) ~ 0(1/L). The energy density 
goes to zero in the infinite volume limit. In this case the 
appearance of the Poincare divergence in Eq. (|80|l has no 
effect on the particle. 

For quantum mechanics the situation is different. We 
can have fluctuations even in non-extensive situations 
|15| due to vacuum effects. For example we obtain, for a 
two-level atom, an energy fluctuation of the dressed ex- 
cited state which is of the order of the decay rate. This 
gives an uncertainty relation between energy and lifetime. 

Coming back to our main discussion, we conclude that 
A^qlqi cannot be expressed as the product Eq. (|79(l since 
A is, by definition, analytic in the coupling constant. To 
make this transformed product analytic, we make the 
replacement 



A- 



E 



c/cl q k Qk 



A 2 ^2 



(84) 



This shows the non-distributive property of A. 

For weak coupling the approximate value of b k is given 
by 0, 



2tt 



L [(w fc -tDi) 2 +7 2 



(90) 



This has a sharp peak at uj k — u>i with a width 7. It 
corresponds to the line shape of emission and absorption 
of the field by the renormalized particle. 

To find more general transformed products hJq\ m q™, 
we apply the same logic that led to Eq. (|86|l . Whenever 
|c;t| 2 appears in Q* m Qi , we replace it with This leads 
to (see Appendix |HJ. 



A+ 



min(m,n) 

E 

a=0 



(91) 



).!«! 



(m — a)\{n — a)!a! 



Q*m-aQn-aya 



where min(m, n) is the smaller of m, n, and 
Y = ^b k q* k q k . 



(92) 



where is a suitable analytic function. Due to the re- 
quirements on A stated in the Introduction this function 
is not quite arbitrary. Indeed, in the integrable case the 
term q k q k would appear in U^q^qi as A 2 ^ fc c k qlq k (see 
Eq. (|47H 1 . In the nonintegrable case c k is extended to 
Cfc or et, and becomes complex. Taking into account the 
requirements (1), (3) and (4) in the Introduction we con- 
clude that a suitable extension of c 2 , to the nonintegrable 
case is the linear superposition 



r 4 



1 



(85) 



where r is a complex constant to be determined. The 
relation r + r* = I is the simplest relation which guaran- 
tees £fc reduces to c| in the integrable case [see also the 
comments below Eq. (|FI3I) ]. 
So we have |H 



where 



A^Jgi = QlQi + ^2 h k1*k1k 

k 



6fc = A 2 |JV 1 |(-|c fc | a +a). 



(86) 



(87) 



As shown in Appendix IF1 using the requirement (2) we 
obtain 

r= expHa/2) Ni = lNi{exp{ _ ia) m 
2 cos(a/2) 

giving a concrete form of A in Eq. I|86|l . By including 
the term b k in Eq. H86|l we have removed the Poincare 
divergence in the product of Gamow modes. As a conse- 
quence, 



Note that b k - 0(1/ L). Hence Y - 0(L°) only if the 
field obeys the thermodynamic limit condition, Eq. |7|l. 
Otherwise Y vanishes as 1/L and A^ becomes distribu- 
tive. Also, when there are no resonances, z\ becomes 
real and both b k and Y vanish. Then At reduces to W 
(see Eq. (jHUJ). In short, both thermodynamic limit and 
resonances are necessary to obtain non distributivity of 
A T in Eq. (|9I|I . which, as we will see in the next section, 
gives the fluctuations found in Brownian motion. 

For A~ 1 q^ m q™ we obtain the expression (|9I|I with 
Qi,Q* replaced by Qi,Q\, respectively. 



B. Obtaining closed Markovian equations 

The A transformation we have presented satisfies all 
our requirements (l)-(4) stated in the Introduction. Now 
we show that A also satisfies the requirement (5) namely, 
that A gives closed Markovian equations. To obtain 
closed Markovian kinetic equations, we first operate A 
on the Liouville equation, to obtain Eq. (JjJJ. Kinetic 
equations involve a projection (or a part) of the ensem- 
ble p. In order for A to give closed kinetic of equations, 
we require that the transformed Liouvillian 9 in Eq. © 
leaves subspaces corresponding to projections of p invari- 
ant. We will represent these subspaces by projection op- 
erators P( v \ which are orthogonal and complete in the 
domain of 9, 



A^Ui ¥= (Atg*)(At 9l ). 



(89) 



The invariance property of 9 is 

p(")0 = 0>("). 



(") 



I. 



(93) 



(94) 
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Thanks to this commutation property, we obtain from 
Eq. (JSJ closed Markovian equations for the projections of 
P, 

i^P {v) p(t) = 9PWp(t). (95) 

We choose pM as eigenprojectors of Lq. We have 
LqP^ — uA^pO) where u>(") are the eigenvalues. Then, 
for the integrable case the relation (|94|l is automatically 
satisfied, since in this case 9 reduces to Lq, the renormal- 
ized free Liouvillian with eigenprojectors p("\ 

In the Friedrichs model the P'"J subspaces consist of 
monomials (or superposition of monomials) of field and 
particle modes. For example the monomials 

q*q k = P {lk) q{q k , tffttfffi = P {lk) qUkQ* qi (96) 

belong to the same subspace p( lfe ) with eigenvalue 
w ( lk ) = C j 1 — 

One may introduce a Hilbert space structure for the 
eigenf unctions of Lq, including suitable normalization 
constants in the Segal-Bargmann representation [2^ ■ We 
also note that transformed variables A^P^A generate 
the kinetic equation J^J as we have ((A^P^ A\p(t))) = 
({A\pM\p)). 

Now we verify that the relation (|94|l is satisfied for 
the A transformation we have constructed. We restrict 
ourselves to the components associated with the mono- 
mials in Eq. (|78fl . These belong to eigenspaces of Lq 
with eigenvalues (m — n)u>i. We denote the correspond- 
ing projectors as p(" m ). Using Eq. (IT4fl in Appendix [i] 
with q[ = 0, we find 

&ql m q™ = l(mz* 1 -nz 1 )q* 1 q 1 -2i 1 mnY}q* 1 m - 1 ql- 1 (97) 

and similarly 

9q* m q^ = [{mz 1 -nz* 1 )q* 1 q 1 + 2i 1 mnY}q* 1 m - 1 q 7 >- 1 . (98) 

Both the l.h.s. and the r.h.s. of these two equations 
belong to the same eigenspace 

p(mn)_ Thif , illustrates 

the statement that 9 leaves the subspaces P'"' invariant, 
satisfying the requirement (5). 

Due to the Y term, q\ m qi are not eigenf unctions of 9, 
so p( mn ) is not an eigenprojector of 9. This is quite nat- 
ural, since the kinetic processes include both the decay 
of the particle modes (through emission of the field) and 
the absorption of field modes. These correspond to the 
first and second terms inside brackets in Eqs. I|97|). (JSSJ, 
respectively. 

C. II subdynamics 

We comment on the theory of II subdynamics, devel- 
oped by the Brussels school We introduce the pro- 
jectors 

n<"' = A^pMA. (99) 



Similar to Eq. i|93[l . they are orthogonal and complete. 
From Eq. £g) it follows that U^L H = L H U^\ Hence 
these projectors define independent subspaces that fol- 
low independent, closed dynamics. The projectors them- 
selves may be written in terms of generalized eigenstates 
of Lh, which give a complex spectral representation of 
this operator [l4| . 

From the completeness relation of the P'"' or IT") pro- 
jectors, we can recover the time evolution in the original 
variables as 

P {t) = ]T A-Wflt) = ]T n^ P (t). (ioo) 

This shows that, as pointed out in the Introduction, dy- 
namics is decomposed into a set of components obeying 
Markovian equations. In order to actually use Eq. I|100|) 
we need to know all the components P^pit) (or equiv- 
alently, U\ u 'p(t)). At the present moment we have only 
obtained a restricted set of these. This is enough for our 
present goal. Rather than solving the original equations 
of motion (which can be done by other standard meth- 
ods) our goal is to show that Brownian motion, Gaussian 
white noise and damping are part of dynamics seen in the 
A representation. 

The contruction of A we have presented here is based 
on renormalized particle modes. A more general con- 
struction of A starts with the commutation relation l|94() 
together with the other requirements stated in the In- 
troduction. The main idea is to associate a "degree of 
correlation" with each subspace p(") . Dynamics induces 
transitions among different P^ subspaces. We have a 
"dynamics of correlations" 35]. This allows us to per- 
form the regularization of denominators of U in a system- 
atic way, depending on types of transitions (from lower 
to higher correlations or vice versa), which leads to A. 
The interested reader can find a presentation of this for- 
mulation in Refs. [H[ll HI. 

VII. COMPARISON WITH THE 
PHENOMENOLOGICAL LANGEVIN EQUATION 

In this Section we discuss the relation between the so- 
lution of the phenomcnological Langevin equation for the 
Friedrichs model and A^ transformed particle modes (i.e. 
dressed modes) . We will focus on the A^ transformation, 
so that the transformed variables decay for t > 01 ( see 
Eq. I|tj9|) ). Remarkably, the Langevin equation and the 
dynamical equations for dressed modes have the same 
solution. 

The phenomcnological Langevin equation for the 
Brownian harmonic oscillator appropiate for the 
Friedrichs model has the form (t > 0) 

j t qx{t) = -iz 1 q 1 (t) + R(t), (101) 

where R(t) is a complex noise source (see Eq. J340 ). We 
use the hats to denote phenomenological variables. The 
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complex coefficient z\ = (b\ — 17 gives the frequency and 
the damping rate of the oscillator. We assume R(t) has 
the following properties: 

(1) R*(t) and R(t') have the delta function (white 
noise) correlation. 



R*{t)R(t')\ = R 2 c 8{t-t') 



(102) 



where ( ) means noise average. 

(2) R(t) has the Gaussian property 



R*(t 1 )...R*(t m )R(t' 1 )...R(0) (103) 

= s mn (R*(tij^t , j j)---(R*(t in )R(t' jn y i 

all pairs 



The noise constant R c is determined in Appendix fus- 
ing equipartition of energy and assuming the bath is at 
temperature T. The result is R 2 — 2 A fksT /Cot. 
Eq. (|101|l corresponds to the equations 

±-x x {t) = -jxi(t) + ^+A(t), (104) 
at m 

jMt) = -^px{t)-rhwlxx{t)+B{t). (105) 

where A(t) and B(t) are independent Gaussian white 
noises [36| . These equations describe a damped harmonic 
oscillator with random momentum and force terms A(t) 
and B(t), respectively. The equations are symmetrical 
under rescaled position and momentum exchange, which 
is consistent with the same symmetry of the Hamiltonian. 
The solution of Eq. HlOlj) is given by 



where 



<Zl(*) = Qla(t) + qir{t), 



q la (t) = gi(Q)e- tfl< 



q lr (t) = er lZlt / dt' R{t')e lZlt 
Jo 



(106) 

(107) 
(108) 



The term qi a (t) describes the damped harmonic oscil- 
lator without noise, and the term qi r (t) describes the 
behavior due to the noise. 

For later comparison we calculate the autocorrelation 
function (q* m {t)q^(t)). We have 

(rw?i(*)) 

= m a (t)+q*i r (t)r(q la (t)+q lr (t)r) 

m n , , 

x - x ml n\ 

= 2^2^ (m-k)\k\ (n-l)\l\ 

k=o 1=0 y ' K ' 

Xfl?- k (t)q?a l (t) (41r(t)ql(t)) • (109) 



The quantity (<7*r(i)<zir(0) is non-zero only when k — 
I, as we can see from Eq. <|103[) . Considering the 



fact that the number of sets of all possible pairs in 
R*(t 1 )...R*{ti)R(t' 1 )...R{t' l )^ is II, we have 

(#*(*)&.(*)) 

x e - tzit / dt 2 R(t 2 y zit2 
Jo 

l 



IIS 



Rt{l-e-^) 



2l 



(110) 



Substituting Eq. (|TU7jl and Eq. ((TTUjl into Eq. l(TTO|) . we 

get 



(qrmut)) 



a i(mz^—nzi )t 



i(m,n) 



E 

1=0 



(m- l)\(n- l)W. 



xg 1 ™^(0)9r i (0)(^)'(e 2 ^ - 1) ; . (Ill) 

Now we can compare the above expression with the 
time-evolved dressed products e !iff *A^q* m g™. We have 
(see Eq. @Hl and Eq. JSTJ) 



e iLBt A 1 q *m q n = 



j.!r). 



(m — a)!(n — a)!a! 



xe i(mzj -nz±)t Q*m-a Qii— aya (112) 



Writing 



=27 



at = E 



1=0 



l\{a-l)\ 



(e 2 ^ - l) 1 



(113) 



and I' = a — I we have 

n n— I 

ALrrt A t „ n 



m\n 



L H t a f *m n _ \ "* \ 

«i » r m _ i _ iyi n _ i _ V \\ m 

1=0 i'=o y ' y ' 

y^gi(mz*—n,zi)tQ*im—l—lQ'n l —l—l'yl-\-l ^27i 

n . . 

mini 



E 

1=0 
n — l 

E 



i(mz* —nz\)t 



{m-l)\{n-l)\l\ 

{m-l)\(n-l)\ 
(m — I — l')\(n — I — ' 



j*m~l — l'Qn—l — l'Yl f 



l'=0 

xY l {e 2 ^ -l) 1 . 
Using Eq. (|!JT|) again we obtain 



(114) 



e iL H t A i q *m q n 



a i(mz^ —7izi)t 



E 

1=0 



7Jr7.l 



(m-l)\(n-l)U\ 



x(A t g* m - i g "- i )r'(e 2 ^ - 1)'. 



(115) 
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Comparing Eq. 1111(1 and Eq. (|115f) . we see the direct 
correspondences 



z\ z\ 



k B T 



(116) 



The form and time evolution of the ensemble average of 
Langevin equation variables are the same as those of A 
transformed variables. Furthermore, if we take the en- 
semble average of A^q™g™, we see a closer correspon- 
dence. Let us assume that the field action J k obeys the 
unperturbed Gibbs distribution. The initial distribution 
p (T) has the form (with (3 = l/(k B T)) 

p Q {T) = Cpoi(Ji,ai)exp(-/3^c; fe J fc ) (117) 

k 

where C is a normalization constant, fcs is Boltzmann's 
constant and T is the temperature. As shown in Eq. (|E5|) 
the average of J k over this ensemble is 



(118) 



To calculate J^k^k (Jk), we need the form of b k . The 
approximate value of b k is given in Eq. (|9U[) . which for the 
weak coupling case is approximated by the delta function 
(2ir/L)5(ujk — u>i) So we get 



^ bk (Jk) = y bJ? B T ' 



k 



k 



(119) 



Note that w^ -1 does not make any divergence for small 
k since bk is proportional to v\ ~ u>k for small k. 

In short, we obtain a complete correspondence be- 
tween A transformed modes and Langevin modes (see 
Eq. Ullfijl ). The systematic removal of Poincare diver- 
gences in the A transformation gives the Gaussian white 
noise structure. 



VIII. THE FOKKER-PLANCK EQUATION 

Using the above results we can now derive the Fokker- 
Planck equation for the transformed density function p = 
Ap. We start with the transformed equation (see Eq. ©) 



d . 

l at p 



(120) 



We derive the Fokker-Planck equation for qi, q\. We 
follow the standard derivation found in textbooks (see 
HUES)- but now in terms of A. Consider a test func- 
tion G(qi,ql), which is smooth and vanishes at \qi \ = oo. 
Multiplying this on both sides of Eq. (|120fl and integrat- 
ing over the phase space, we have 



d 

dTGiq^qDi— p(T,t) 



dTG(qi,qt)e(T)p(T,t) 



= / drdr'G(q 1 ,qi)0(T)6(T-T')p(r',t). 



(121) 



In Eq. i|121|) . 0(T) means 9 that acts on T variables. We 
expand G(qi,ql) near q[ and q{* . 



J dTG{q x ,ql)i^- t ~p{T,t) 



(122) 



= I dTdT'l Y y n ° °,. G(q!,g>) 



x(«I-? 1 *r(?i-^) n jfl(W-r / )p(r / ,*). 

Integrating by parts, Eq. I|122|) becomes 

j dTG(q u ql)i^p(r,t) (123) 



= / dT'G( q [ iqi *)J2f2 



^ ^ m\n\ d(q'{) m d(q[y 

m-0 n-0 v 1 J v A/ 



dT(q* 1 -q 1 *) m (q 1 -q' i re(T)S(r-r') 



p(T',t). 



We call the quantities inside the brackets in Eq. (|123[1 the 
"moments" of order m + n. The moments are calculated 
explicitly in Appendix [I] They are given by 



dr(ql- qi *r(qi-q[) n mS^^^) 



ziq[, m = 0, n = 1 
— z\q*, m = 1, n = 

2i lJ2k b kQk<lk^ m = 1,71 = 1 
0, for all other m and n. 

Substituting Eq. (1 12411 into Eq. I|123|l . we get 

j dVG(q 1 ,q* 1 )z^- t p{T 1 t) 



(124) 



(125) 



dT'G(q' 1 ,q 1 *)(- — )(z 1 q' 1 )p(T',t) 
dq[ 

+ I dr'G(q' 1 ,q 1 *)(-^)(-zfqnP^\t) 

oqidqC 

By changing the integration variable T' to T in the 
right hand side of Eq. I|125|) and eliminating i on both 
sides, we have 

J dTG{ qi ,ql)j t p{T,t) 



dqidql 



(126) 



Now suppose that p(T,t) is factorized at t — 0. In other 
words, we write p(T, 0) as 



p(r,0) = gi{qi,qi)Y[gk{qk,q* k ) 



(127) 
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As shown in Appendix|jJ this factorized form of p enables 
us to write Eq. H126|) as 



J dTiGfo.gJ)^ J dT f p{T,t) 



dr 1 G(q 1 ,q* 1 ){^-(iz 1 q 1 ) + jL(-iz*q*) 
oqi dq 1 

— (2 7 J2 b k(<lUk)) J dT f p(T,t). (128) 



Since G(ri) is an arbitrary test function, we can write 
Eq. QTM as 



m hFi,t) = {izr-^-qr-izt—qt (129) 



where 



pi{V u t) = J dT f p(T,t). 



(130) 



Eq. (|129l) is our Fokker-Planck equation for the normal 
modes. This equation is applicable for any initial field 
configuration obeying the extensive condition and Eq. 
(|127fl . In the non-extensive case, the diffusion term con- 
taining bk vanishes, and the equation describes damp- 
ing of the oscillator without Brownian motion. For the 
special case where the field has the unperturbed Gibbs 
distribution, using the approximation (|119fl we get 



2 1 k B T d 2 
U>1 dqidql 



}p 1 (r 1 ,t). (i3i) 



This is precisely the equation one obtains from the phe- 
nomenological Langevin equation 1)101(1 . 

The Fokker-Planck equation for other variables can be 
also derived from Eq. (|129fl by changing variables. For 
example, the Fokker-Planck equation for the position and 
momentum x\ and p\ is given by 

— pi{Ti,t) = \ - — (— - 7x1) + —(mufx-L +7Pi) 
at L oxi m opi 

D B 2 D 3 2 -1 



where 



fa = muji/uii, 

d x = j^- y~] bk (Jk) ~ ^ lkB J , 

mu>i *■ — ' mujt 

k 

Dp — 2ifi"fCji bk (Jk) ~ 2rh^fkBT. (133) 



[the approximate values are applicable for the unper- 
turbed Gibbs distribution.] The Fokker-Planck equation 
for the action variable J\ is given (after integration over 
the angle variable u\) by 



(134) 



D;, 



d 2 
dJ 2 



where 



Di 



k 



2jk B T 
&1 



(135) 



Eqs. I|132|) and (|134|l coincide (in the weak-coupling 
approximation) with the equations for Brownian motion 
of an oscillator in an anharmonic lattice derived in Ref. 
35] . Eq. 11341) (in its exact form) was first proposed by 
T. Petrosky J33f- 

Note that Eq. (|132|l is symmetric with respect to 
rescaled position x\ and momentum p\. The reason is 
that the Hamiltonian considered here is symmetric in 
rescaled x\ and p\ to begin with. The same is true for 
the anharmonic lattice model considered in Ref. |3^. In 
contrast, the Kramers (Fokker-Planck) equation |35l l38l ] 
derived from the Ornstein-Uhlcnbcck phcnomcnological 
theory of Brownian motion [39| is not symmetric, because 
the Brownian force breaks the position-momentum sym- 
metry. In spite of the difference, for the case 7 <C uj%, Eq. 
(|132fl gives the same solution as the Kramers e qua tion 
The solutions of Eq. i|132|) can be found in Ref. 



IX. CONCLUSIONS 



quat 
35]. 



In this paper we studied the irreversible and stochastic 
behavior of an oscillator coupled to a field in the ther- 
modynamic limit, using the star- unitary transformation 
A. We showed that the average of dressed particle modes 
has the same time evolution as the ensemble average of 
particle modes in the Langevin equation (|101|) . Also, 
the reduced distribution function for the particle vari- 
ables exactly obeys the Fokker-Planck equation (|129fl . 
which describes the damping and diffusion processes. It 
is remarkable that the systematic removal of Poincare 
divergences by analytic continuation leads to the same 
structure as that of Gaussian white noise. 

Since the Gaussian structure of the fluctuations is com- 
ing from the resonances, rather than the specific form of 
the initial field ensemble, our derivation of the Fokker- 
Planck equation is valid for both Gaussian and non- 
Gaussian field ensembles. Due to the Fokker-Planck 
evolution, the particle distribution is still Gaussian for 
t — > 00. 

Our method of isolating Poincare divergences can in 
principle be applied to more complicated systems than 
the one considered here, in order to investigate further 
the relation between noise and dynamics. 
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We have studied the fluctuations of a particle sur- 
rounded by a field. One can also consider the fluctuations 
of the field modes induced by the particle. An interest- 
ing result is that fluctuations of the A-transformed field 
modes appear even when the energy of the field is non- 
extensive. If we go back to the initial formulation, acting 
with A -1 , we recover the usual equations. This will be 
discussed elsewhere p3 |. 

Markovian equations such as the Langevin or Fokker- 
Planck equations describe irreversible and stochastic pro- 
cesses. Rather than viewing these equations as approx- 
imations of Hamiltonian dynamics, we see them as de- 
scribing components of dynamics. These components can 
be identified through a change of phase-space variables 
obtained by the A transformation. This gives a represen- 
tation with broken time-symmetry. 

In the non-Markovian equations, the effects of dressing 
are not separated from irreversible (or thermodynamic) 
processes. In our approach the dressing on the particle 
is incorporated from the beginning, since we deal with 
renormalized particles. This allows us to isolate pure 
thermodynamic behavior such as Brownian motion with 
a white noise source. In this sense, we can aim to refor- 
mulate thermodynamics in terms of renormalizedparti- 
cles. This is interesting in view of recent claims |j| that 
traditional thermodynamics may not apply for quantum 
Brownian motion at low temperatures, due to the non- 
Markovian character of quantum noise. 

The relation with thermodynamics requires a precise 
definition of entropy or 7i function. The latter can be 
given in terms of A ^(|. The introduction of A leads 
to probabilistic considerations in classical mechanics, in- 
dependent of incomplete knowledge or quantum correc- 
tions. In conclusion, we believe that our approach leads 
to a unification of dynamics, thermodynamics and noise. 



APPENDIX A: EXACT TIME EVOLUTION OF 
PARTICLE AND FIELD MODES 

We consider first the integrable case. The time evolu- 
tion of the modes q s (t) can be calculated using the renor- 
malized modes Q s , which are eigenstates of the Liouville 
operator. We have 



(Al) 



The particle mode Q\ is given Eq. (|47l) . The field modes 
Qk are found from the equation LhQr = —&kQk, which 
gives 



Qk = K 

where 



1/2 



qk 



AVfc 



fe+ E 



_XVk_ 

U>k — LOk' 



r) k (z) = z - uji - } t 



fe'(^fe) 



z - uj k > 



qw) (A2) 



(A3) 



ui k = ui k + 



(A4) 



N k = 



Vli^k) 



d+ E 



A 2 ^ 



k'(^k) 



(uj k ~ u k >f 



(A5) 



Note that Co k = u k + 0(1/ L) and also N k = 1 + 0(1/ L). 
As long as L is finite all the denominators are non- 
vanishing, and there are no resonances. 

We write the linear relations between dressed and bare 
modes as 



Qs — ^ ^ ^-sr^ri Qs ^ ^ CrsQr 



(A6) 
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Then, using Eq. (|A1(I we get the coefficients f sr (t) in the 
equation q s (t) = Y, r fsr(t)qr(ty as 



fsrft) ^ c rr i e r c r / s 



(A7) 



In the nonintegrable case resonances appear (see Sec. 
Ivjl and if we insist on keeping the renormalized modes as 
usual canonical variables, then the particle modes dissap- 
pear into the continuum of field modes (we can however 
introduce the Gamow modes of Sec. which are non- 
canonical). Keepin g ca nonical modes, the Hamiltonian 
is represented as 0, [lj, l22l | 



H = J2^kQ* k Qk 



(A8) 



in the continuous limit. The renormalized field modes 
are given by 



Qk 



AU, 



qk 



V T (uk) 



Qi+ V —qk' 



(A9) 
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where e is a positive infinitesimal quantity. There are two 
branches, corresponding to analytical continuation to the 
lower or upper half planes of u> k . For t > 0, and given 
the initial condition Eq. (|41|) , it is convenient to take the 
branch with — ie, since this will give exponential decay of 
the particle modes the positive t direction. We get 



(*) - E -S^-Mv = E -#rft(^" lM 



- I 4. 1 v l2 -9i(°K 



EA 2 VJ 2 x ^ 



, ,?7 + (^fc)h '-g' - Uk> - it 



q k >(0)e 



-iuj k t 



(A10) 



For t < 0, we take the branch with +ie. 



APPENDIX B: RANDOMNESS IN THE FIELD 
MODES 

We choose our initial condition with the form (|41|l . 
In classical mechanics qio = (?i(0) can be determined 
exactly since q±o is a function of the initial position and 
momentum of the particle. For the modes q k Q we need 
more care. Let us first write q k Q in terms of action and 
angle variables, 



q k o = V^e~ Mfc °- 
With the ensemble 141(1 we have 



lim (Jfco) 

L — >oc 



O(L ) 



(Bl) 



(B2) 



in the thermodynamic limit. For example, for an unper- 
turbed Gibbs thermal distribution of the field modes we 
have (J k0 ) = k B T/uj k (sec Eq. (|E5|) ~). 

For almost all phase points { J10, --Jko--, ■•, •■} 
out of the ensemble, any two different angles a k o an d 
a k >o have no correlation. In other words, the sequence of 
angles {ce kn o} is completely random for almost all cases. 
This property allows the existence of the thermodynamic 
limit [l8j. Indeed, if a k ^ is uniformly distributed over 
[—it, 7r] and the sequence of angles {ckfc„o} is completely 
random, then the term 



E 



E 



AVfc 



/ 7, „-i(^k+ak(\)t 
J k Q c 



(B3) 



in Eq. (|A10|I is O(L ). This is so because the summation 
is taken over a sequence complex numbers with random 
phases and hence it is proportional to the square root 
of the number of modes, which in turn is proportional 
to L. Since V k ~ L" 1 / 2 , Eq. (£52) is 0{L°). It shows 
a very irregular time evolution as the number of modes 
increases. 



Note that if a k o was a smooth function of k, then for 
the first term of Eq. i|ATfJ|) we would have in the limit 
L — > 00 

XV k 

I . — T7 ^QkO 

V V + {uk) 

L 2ir m—^ Xv k 
9Z~T 2^' 



2tt L r)+(cj k ) 



J k oe 



s/*5TO^'"*" (B4) 



and since the integral is O(l), this expression would di- 
verge as 0(yL). 



APPENDIX C: DERIVATION OF THE 
NON-MARKOVIAN FOKKER-PLANCK 
EQUATION 



From the Liouville equation we have 



J dTG(T 1 )^- t p(T,t)=-i J dTG(Ti)LHp(T, t) 

OO OO 

= EE Gm « M ™«w ( C1 ) 



m— n— 



where 



M mn (t) = -if dTql m q™L H p{T, t). (C2) 
Using Eq. (|5U(l we have 

M mn (t) = { dT{iL H [ql{t)] m [ qi {t)] n )p{T,Q). (C3) 



Since Lh is a differential operator, we have 

iL H [qt(t)r[ qi (t)} n = m[?J(*)] tn - 1 [gi(t)] n ti ff ?J(*) 
+ n[qi(t)] m [ qi (t)] n -HL Hqi (t). 

(C4) 

From Eqs. 1(31(1 and ((34(1 we have 

iL H ql{t)=iz$(t)ql(t)+R*(t), 

iL Hqi (t) = -i Zl (t) qi {t) + R(t). (C5) 

Inserting Eqs. (fCljl and lfC5)l in Eq. l(C3jl we get 

M mn {t)=M* mn {t) + M* n {t) (C6) 

where 

M* n (t) (C7) 
<ZT[«; (t)m - i Zl (t)n\ [ql (t)) m [ qi (t)] n p(T, 0) 



= / <T[iz%(t)m-izi(t)n]q? n tfpQ? ) t) 
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and 



M* n (t) = M^(t) + M%*(t), (C8) 
= [ drm[ql{t)] m -\ qi {t)] n R*{t)p{T,Q), 



M^(t)= J drn[ql(t)) m [ qi (t)} n - l R(t)p(T,0). 

Now we evaluate M%£(t) in Eq. fHsjl . Using Eq. we 
have 

(C9) 

dr mfeKi)]™- 1 [<&(*)]" E ^(*M/>(r, °)- 

Then, using Eq. we obtain 



M ' (t) 



dT 771 



/u(*)9i + X)/ip'(*) 



9 P ' 



p' 



(CIO) 



5>fe(*Kp(r,o). 



Abbreviating / Q = fi a (t), hk(t) = /ife and using binomial 
expansions, we get 



m — 1 n 



(m — 1)! n! 



_ n ; Q (m-l-a)lo! (»-6)!M 
xm[/rg 1 *] m - 1 - Q [/ 1 g 1 ]"- b (Cll) 

X E E fpi" ' fpa^Pa+lfp'l" ' fp' t 
PX-Pa+1 p' x ~p' h 

X 1pi ■ ■ ■ Ipalpa + llp'l "'IP', ^( r ' ) 

where we changed the variable fc to p a +i Due to the as- 
sumed form of the field distribution in Eq. (|41|l . only 
expectation values of observables independent of the an- 
gles a p of the field modes are nonzero (in other words, 
each q p > field mode must be multiplied by its complex 
conjugate q*. with p'j = pi). Thus in Eq. (|C11(1 we must 
have b = a + 1 and we get (a + 1)! possible parings 

(m— X,n— 1) 



M^(t) = [dT £ 
(to - 1)! 



(m - 1 - o)!a! (n - 1 - a)!(a + 1)! 

xTO[/r <z *] m - i - a [/ igi ]"- i - a (ci2) 
x E i/ P1 i 2 ---i/ P JX +1 /^ +1 

x(a + l)!g* i g Pl • ■■q* Pa q Pa q* Pa+1 q* Pa+1 p{T 1 Q). 

Note that due to the volume dependence of the coeffi- 
cients f p ~ 0(1/Vlj) and /i p ~ 1/VL [which follows 



from Eq. (|AT0)| ]. we can neglect the cases where two or 
more of the wave numbers are repeated, e.g., pi = Pj for 
i j. Indeed, these cases give contributions of order 
0(1/L) within the summations in Eq. (|C12() . Canceling 
(a + 1)! and writing k — p a +i and n\ = n(n — 1) we get 



M J 



(C13) 



/mm(m-l,n-l) 
dr e 

x (m-1)! (n-1)! 
(to — 1 — a)!a! (n — 1 — a)! 

x [/rg 1 i m - l - a [/i9i]"- l - a 
x E \u 2 ---\f P J 2 

Pl—Pa 

XQ^Qpi ■ ■ ■ Qp a Q Pa E h *kfk<lUkp(T, 0). 



Applying the reverse steps from Eqs. (|C12|) to HC10|) we 

get 

M^(t) = ran J dT [tf (t)]"- 1 ^)]*- 1 (C14) 
x E^(*)/ifc(*)9fc9fcp(r,0) 

= to 7i | dTKCt)]™- 1 [gi(f)]"- 1 

x $>iEW/i*(*) («**) p(r,o) + o(i/L) 

k 

where due to the factorization property Eq. 141|) we can 
take independently the average 

(€%) = (Jk) = J dTqtq kP {T,Q). (C15) 

[a similar argument is given in Appendix \j\ . Bringing 
the time dependence back to p we get then 

x ^h%(t)fi k (t) (J k ) P (T,t). (C16) 



For the second term in Eq. (|C8I) we have M^(t) = 
[M^[t)\* . Putting everything together in Eq. IjCip with 
integration by parts and a few straightforward manipu- 
lations, we obtain 

J dTGiT^p^t) 
= J dTG(Ti) " <*i(*)^T?i 



D(t) 



dqidql 



p(T,t). 



where 



(C17) 



(C18) 
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Since G(Ti) is an arbitrary test function, we come to the 
non-Markovian equation i|42|) for the reduced distribu- 
tion function pi. Replacing the explicit form of h k (t) we 
obtain as well Eq. (|43|l . 



APPENDIX D: MARKOVIAN APPROXIMATION 

Phenomenological equations may be obtained from dy- 
namics using the so-called Markovian approximation (or 
the X 2 t approx imation studied by Van Hove and others 
[37L l40Ll4ll|42ll) where memory effects are neglected (see 
also [43I I44L l45j p. This approximation is valid if we take 
the continuous spectrum limit L — > 00 (with Poincare 
resonances), for weak coupling and for time scales where 
the relaxation process dominates over dressing pro cesses 
(the particle-bath correlations are neglig ible |4|,|35|). 

In this appendix we consider the Markovian approx- 
imation of the dynamical evolution of the bare particle 
modes. We will write the weak coupling approximation 
of the coefficients fu(t) and fik(t), which will give the 
Markovian approximation of Eqs. (|34|) and 142|) . 

We start with Eq. I|A10|) . for t > 0. We approximate 
for weak coupling A< 1, 



1 



1 



r/+(uj k ) wu-Zi 



(Dl) 



We then separate the pole contribution at ujk = z\ 
and the branch cut contribution from each term. The 
pole contribution gives the exponential decaying part and 
the cut contribution give s classical Zeno effect and non- 
exponential behavior [22j • Here we will only consider the 
pole contributions in Eq. IjAlOfl . which amounts to the 
Markovian approximation. Taking the pole contribution 
at uik — zi in the last two terms of Eq. i|A10[l , we get 



qi (t) « qi (o) e - i2it + ]T 
_y^_^_ gfc(Q)e -i,xt < 



AVfc 

UJk - Zl 



■9fc(0)e 



— lUJ^t 



(D2) 



Therefore we have 

/n(*)«e 



— iz\t 



hk{t) « -^Me-*"** - e" 4 "- 1 '! 



(D3) 



This leads to 



zi(t) 
R(t): 



E- 



/ > ' XVke-^qkiO). 



(D4) 



For t < we can repeat the same procedure, except that 
now we choose the +ie branch in Eq. I)A9|I . Then we 
obtain the following approximate equation 



qi(t) w -iz x qi(t) +R(t), t>0 
qi(t) w -iz*qi{t) + R(t), t<0 



(D5) 



which has a Langevin form. This equation is expressed in 
two branches, one for t > and the other for t < 0. The 
change —izi => —iz* corresponds to the change 7 =>■ —7, 
since Zi = u)i — 17. Taken separately, each branch breaks 
time reversal invariance, while as a whole time reversal 
invariance is kept [note that in the integrable case we 
have 7 = and there is no splitting into two branches]. 
Comparing Eq. (|D5(1 with Eq. (|101fl we identify zi = zi . 
Furthermore, the term R(t) behaves as the white-noise 
source R(t) in the sense that in the pole approximation, 
the auto-correlation functions of the variables qi and qi 
coincide. Indeed, the "noise" term R(t) in Eq. 1D4I) has 
the same Gaussian property (| 1 03|) as the noise R(t), pro- 
vided we replace noise averages with averages with en- 
sembles of the form ()41|l (see discussion below Eq. i(CTT|l). 
Then, defining 



qi a (t) = qi(0)e- 1 ^ 



(D6) 



qir(t) = e 



dt'R(t')e l 



AV fc 



uj k - zi 



(D7) 



and following the same steps as in Eqs. (|106fl - l|110|) we 

get 



(qt;(t)q[ r {t)) = USki (qUthiAt)} 1 
In general we have (with q k Q = <7fe(0)) 
(qt(t + r)q lr (t)) 



(D8) 



E 



uj k - z 



7 * ( P iu k (t+T) _ e izl(t+T)\ 



E 



UJl - Zl 



(D9) 



For the normalized thermal field ensemble we have 

k B T 



(QkoQk'o) = Skk 1 



UJk 



(D10) 



Using this result and going to the continuous limit we 
obtain from Eq. l|D9ll 



(q* lr (t + T)qi r (t)) 



dw- 



/ \w-zi\* uj 

(Dll) 

For 7 <^ uji the integrand is sharply peaked around uj — 
uji. We separate the pole and the cut contributions to 
the integral, rewriting Eq. HD11|) as 

(q* lr (t + r)qi r (t)) (D12) 

~ (Qlrit + T )<llr(t)) pole + (9lr(* + T )^r{t)) cut 
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where 



{<hr{t + T )Qlr(t)) pole 



\u> — Z\ | 2 U! 

(D13) 



(q*lr(t + r)qir(t)) cut 



\ 2 vl k B T 



\ui — Zi\ 2 UJ 

( e iuJT — e -izit e iu(t+r) _ e iz'*(t+T) e -iut _|_ e iz i T g- 2 7t^ 

(D14) 



Using 



Vzi ~ VQ lt 



k B T k B T 



7 « ^A 2 ^, (D15) 
the pole contributions at u> = zi, u> = z* give 

(q$ r (t + r)q lr (t)) pole « ^>i T (l - e" 2 *). (D16) 

The cut part involves the tail of the Lorentzian distribu- 
tion l/\w — zi\ 2 — 1/({lu — uoi) 2 + 7 2 ). It gives 

{q{ r {t + T)q lr (t)) cut 

X 2 v 2 k R T 



dui- 



l \U3 + Zl\* lu 

x (eT i ^ r — e - iz i t e - iuJ ( t + T ) — e iz i(*+ r ) e iwt _|_ e iz i T e~ 2jt ) 

(D17) 

In weak coupling case (7 <§; the cut contribution is 
much smaller than the pole contribution. Replacing the 
result (|D16(I in Eq. (|DB|> with r = we obtain the same 
correlation Eq. (|lllf) obtained from the phenomenologi- 
cal Langevin equation. 
We have also 



D(t) 



2^k B T 



(D18) 



Hence we recover the Markovian Fokker-Planck equation 
(|131|l with Z\ = Z\ . 



APPENDIX E: PROOF OF EQ. dsTl) 

In this appendix we show, by perturbation expansion 
that Eq. 181|) is satisfied for the equilibrium Gibbs en- 
semble. The same relation is valid for the class of ensem- 
bles (generally non-equilibrium ensembles) having delta- 
singularities in the wave numbers, of which the Gibbs 



ensemble is a special case |13l 114 135| . The Gibbs ensem- 
ble is given by 

p cq = Cexp(-(3H) (El) 
where C is a normalization constant, such that 

J dTp cq (T) = 1. (E2) 

In the perturbation expansion we have 

p c « = Ccxp(-I3H ) ^1 + XV + ^X 2 V 2 + ■ ■ . (E3) 

Due to the angle integrations, only diagonal monomials 
are nonzero 



JdrH q* r m "q^ exp(-pH Q ) oc [] s. n 



(E4) 



Considering the explicit form of the potential V we then 
have 



(J k ) = J drj kP c * 



J dTJ k ex V {-(3H ) , „ 2 , 



J dTexp(-f3H ) 



0(X Z 



J °° dJ k Jk exp(-/3uj k J k ) ( 



/ °° dJ k exp(-/3w fe J k ) 



x 2 v k v k , 



fdrgl qk ,(X 2 V 2 /2)cxp(-[3H a ) 
J dTexp(-f3H ) 
1 1 1 



f3u>i fiu> k f3u> k > 



0(X 4 



(E5) 

0(A 4 ) 
(E6) 



Noting that V k ~ L 1 / 2 , and Y^y ~ 0(L), we obtain the 
result JEIJ. One can check as well that terms of higher 
order in XV retain the same volume dependences. 

For an ensemble of the form Eq. (|41|l . we have (q^qk') = 
0, while (Jfe) ~ 0(1). This is a special case of Eq. (IFT|) . 



APPENDIX F: A AND PRESERVATION OF THE 
MEASURE 

From the requirement (2) in the Introduction, A pre- 
serves the measure of phase space. This means that 



dTAp = 1, 
dTA^p = 1 



(Fl) 
(F2) 



for any normalized ensemble p. The first equality is easily 
shown, since it may be written as 



J dT{Kh)p= 1 



(F3) 
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where 1 is the unit operator. Since At can be expressed 
as a perturbation expansion, 



A f = 1 + 0(XL V ) 



(F4) 



and Lyl — we conclude that A T 1 = 1, from which Eq. 
(IF1|) follows. Similarly, one can show the second equality 
(lF2l) . 

So, Eq. 1F2J1 should be satisfied for the specific A we 
have constructed in Sec. IVII This condition will allow us 
to derive the relation (|88|l . Consider the ensemble 



p = Ciq{qi exp(- J/ J ) 
where Ci is the normalization factor given by 



C s 



with 



dTq*q s exp(-J/J ) 



(F5) 



(F6) 



(F7) 



and Jo a constant that makes the argument of the ex- 
ponential dimensionless. The factor exp(— J/Jq) en- 
sures the existence of a finite norm of p, (see the Segal- 
Bargmann representation in |23,l49)- The total action J 
is an invariant of motion, because we have Lq J = and 
L V J = 0. Using Eq. (F4j), we get 

At J = J, At cxp(- J/ J ) = exp(- J I J ). (F8) 

The operator Ly is a differential operator. Applying the 
chain rule of differentiation and Eq. (|F8|) we conclude 
that 

AUlqi cxp(- J/ Jo) = (A+glffi) exp(- J/ J ). (F9) 



Inserting the ensemble Eq. (|F5|) in Eq. (|F2|) and using 
Eq. JHEJ) we get 

d f dT\N x \[q{ qi + A 2 ^(r4 + c.c.)q* k q k ] 

^ k 

x exp(- J/ J ) = 1. (F10) 

where the off-diagonal terms such as q^qk appearing in 
the product QlQi in Eq. I|86|l vanish due to the integra- 
tion over angles in phase space. We can write Eq. (|F10|> 
as 

dl^lpf 1 + A 2 ]T(rc 2 + c.c.)^ 1 ] = 1. (Fll) 

k 

Since d = C k for any k (see Eq. (jHl), Eq. (jF10)l leads 



to 



|iV 1 |[l + A 2 E(r C 2 +c.c.)] = l. 



(F12) 



This equation plus the condition r + r* = 1 yield the 
result H88fl . With this result we can write 

A f gjgi = Q (0) Q\Qi + P (0 \rQlQi + c.c.) (F13) 



where is the projector to angle-independent mono- 
mials (such as q^qs) and = 1 — p( '. This shows 
that the transformed product of particle modes can be 
simply written as a combination of renormalized particle 
modes (i.e., Gamow modes). 

The derivation followed here is similar to the derivation 
followed in Ref . , where we used the A transformation 
to define dressed unstable states in quantum mechanics. 
The only difference is that in the relation r + r* = 1 
was derived from the requirement that the dressed un- 
stable state has an energy fluctuation of the order of the 
inverse lifetime. This fluctuation is a quantum effect. 
Here we are dealing with classical mechanics, so we pos- 
tulate r + r* = 1 as a basic condition. An alternative 
derivation, presented in Appendix A of Ref. |lfi| | , started 
with the analogue of Eq. I|F13(1 . as a postulate. All the 
derivations give the same result iJHHJ- Note that a con- 
dition different from r + r* = 1 would not allow us to 
express Atq^ as a combination of Gamow modes only, 
and it would lead to energy fluctuations different from 
the inverse lifetime in the quantum case, which would be 
unphysical. 

The preservation of the measure can be proven for 
more general ensembles involving monomials of particle 
modes, which we have considered in Sec. IVII Indeed, 
from the relation ^2 k c* k Ck = —1 (see 0|) we have (see 
Eqs. (1FT21) . JH3) 



5> = i. 



(F14) 



Using this relation as well as the expression Ij91|l one can 
show that for p = C mn qi m qi exp(— Jj Jo) with C mn a 
normalization constant, we have 



J dTk^p = j dTp = 5„ 



(F15) 



The proof (which we will omit here) uses the relation 
/ jrQf n Q?exp(-J/J ) = 0. This follows from the fact 
that both Ql m Qi and exp(— J/J ) are eigenfunctions of 
Lh with different eigenvalues, which implies their orthog- 
onality. 



APPENDIX G: PROOF OF ANALYTICITY OF A 



In this appendix we show that Eq. (|91|l removes all the 
non-analytic \c k \ 2 terms, replacing them by t; k = rc k +c.c. 
First we derive recursive formulas to calculate A'q* m q". 
We start with Eq. (|9"T|l for m> n 



At 



7.!t7.! 



Q*™-lQn-l Y l. (QI) 



(m-iy.(n-l)W 



(the n > m case can be calculated by taking the complex 
conjugate of Eq. HG1|) 1. We have as well 

Atq* m+1 ^ 
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E 



(to + l)\n\ 



QJ m + 1 -'Q»-'y'.(G2) 



^ (m + l-0!(n-/)Wl 
For / > we have the identity 
(m + 1)! to! 



,(G3) 



(m + l-Z)!Z! (m-0!i! {m + 1 - - 1)1 
Inserting this in Eq. I{G2|I we get 



,1 



TO 



^ L(m-i)W (m+l-0!(i-l)! 

Til 

/o*m+l— lr\n— l\rl 



Q™~ l Y l 



(G4) 



The first term plus the second term give 

(A+?r«i)A t «I (G5) 

(note that A^g* = QJ). The third term may be written 
as (with V = I — 1) 



E 

Z'=0 



to! n(n — 1)! 

(to -I'Wy. (n-l'-iy.^ 1 



Therefore 



A+ ?1 m+1 9i 

= (A^r'^DA^i* + nYA^ql" 1 ^- 1 (to > n) 
For to > n we have, from Eq. (|G1(1 . 



(G6) 
(G7) 



AUl m Qi +1 = E 



m!(n + 1)! 



(m-iy.{n + i - iy.n 

x Ql m - L Q 7 > +1 - l Y l . (G8) 
Using Eq. (|G3[> we get 



A t g *m ? „+1 

= QTQi +1 



(m — n — 1)! 



77i —n—1 y^n-\- 1 



E 



(m-iy. i(n~iyn (n+i-iy.(i-iy. 

xQl m - l Q^ +1 - l Y l . (G9) 
Adding the first and the third terms we get 

(Atgr#)At 9x (G10) 

(note that A^gi = Q\). Adding the second and fourth 
terms we get (with V = I — 1) 



E 



to(to — 1)! 



l'=0 



(m- /'-!)! (n-l')l(l'y^ 1 



*m — l' — lQn — l +1 



toYA 1 ^*™- 1 ^ 



(Gil) 



Therefore 

= (A^^q^Atqi + mFAtg*" 1 -^™ (to > n). 

Eqs. I|G7|) and (|G12I) plus their complex conjugates per- 
mit one to construct Atg" n q" recursively. 

Now we prove the analyticity of A^ql m q" at A = 
from the recursive relations. In the recursive relation, 
we show that if the lower order terms in to and n 
like AU* m qf, Atg™^ 1 - 1 and A^q*™- 1 ^ 1 are analytic, 
then the higher order terms A^* m+1 g™ and A^q^q™ are 
also analytic. Then from mathematical induction, the 
analyticity of A^"™^™ is proved for general to and n 
(the to < n case can be shown in the same way). In Eq. 

13, 



A i q *m+l q n = (A t g 



*m n 



E 

1=0 



m'.n 



(m-iy.(n-l)U\ 



)A t g* + nYA^q^q™- 1 



+nY 



■ n-l 

E ( TO _ >)|( n 



m!(n - 1)! 



1-/)!/ 



\*m — If^n— 1 — lyl 



■QT QT 



(G13) 



Suppose that the quantities inside large parenthesis are 
analytic in A. The additional non-analytic terms appear 
whenever additional products QlQi appear. Since 



Ql = N 1 1 /2 (q 1 



k 



(G14) 



each QlQi produces a |c/c| 2 term, which is non-analytic 
in A. Let us denote the non-analytic part of a function 
/(A) as Fn(f(X)). The non-analytic part in the first 
term in the right hand side of Eq. (|G13|) is made by the 
additional Ql multiplied by which generates n — I 

terms |cj;| 2 : 



Fn 



i.lnl 



(m-iy.(n-l)U\ 



Emlnl 
; _ o (m - iy.(n - iy.n 



Q* m -lQn-l-lyl 



X(n-/)A 2 |AM Y\c k \ 2 q* k q k 



n-l 

E 

1=0 



k 



mini 



Q *m-L Q n-l-lyl 



(m - iy.(n — i — iy.n 



xA 2 |^|^ 



Ck\ 2 qiQk- 



(G15) 



The non-analytic part in the second term in the right 
hand side of Eq. i|G13|) is coming from Y. Since 

Y = E bk 1kQk 

k 

= ^A 2 |iV 1 |(-| Cfc | 2 +rct+r*c^)qiq k , (G16) 
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the non-analytic function \ck\ 2 appears inside Y. 



Fn 



n-l 



nY^ 



m\(n — 1)! 



1=0 



(m-J)!(n-l-0!J! 



m — l^n—l — l^l 



= -#\n x \y, 



E 



1=0 



(m-l)\{n-l-l)\a\ 



Q* m -lQ*n-l-l Y l, (G17) 



The non-analytic parts from the first term and second 
term in Eq. ljG13(l exactly cancels out. So, the left hand 
side of Eq. (|G13J| is analytic in A. Note that terms of the 
form \ck\ 2n with n > 1 give 0(1/L) contributions and 
thus they are negligible. 

Next, we show that the left hand side of Eq. 1|G12|1 is 
analytic in A. The non-analytic part of the first term in 
the right hand side of Eq. (|G12I) is 



Fn 



= Fn 



(Atgr«?)A+5i 

n 

E 



(m- l)l(n-l)U\ 



Q*m-lQn-l Y l)Q 1 



E 

1=0 



{m-l)\(n-l)\l\ 



Q*rn-l-lQn-lyl 



:(m~l)X 2 \N 1 \J2\ck\ 2 q* k q k 

n 

= E 



k 

mini 



1=0 



(m-l- l)!(n-Z)W 



■Ill —l — l^n — lyrl 



a 2 \ Nl \Y, 



Ck\ 2 q* k qk 



(G18) 



The non-analytic part of the second term in the right 
hand side of Eq. 1|G12|) is 



Fn 



(G19) 



E 



1=0 



(m- l)!n! 
(m - I - l)\(n - 



QT 



i-l-?Qn-/yl 



E 



(m - I - l)\(n - l)W 



*m— 1— I s\n— 



Q^~ l Y l 



Again, the non-analytic parts of the first and second 
terms of Eq. I|G12(1 exactly cancel out. The right hand 
side of Eq. (|G12(1 is analytic in A. Therefore from the 
mathematical induction A'g™^ is analytic in A. 



APPENDIX H: CALCULATION OF THE NOISE 
CONSTANT R c 

In this appendix we determine the noise constant R c . 
We assume that the noise R(t) comes from the thermal 
bath with temperature T. In this case, we expect that 
the system reaches thermal equilibrium for t — > oo. Fur- 
thermore, from the equipartition theorem we expect that 



(Pi) eq _ 1 



e i 2m 



knT, 



(HI) 



where kg is Boltzmann's constant. Substituting the re- 
lations 



1 



2m p tl>i 



(&(<)+£(*)), 



m v LO\ 



(fc(t)-«(t)) (H2) 



into Eq. I|H1J| . we get the conditions 

{qi(t)) eQ + mt)) eQ = 0, 
) eq 



eq \ J 1 v ; l eq 

&i (q*i{t)qi(t)) ea = k B T. 



(H3) 
(H4) 



On the other hand we have 



= lim((q* la (t)qia(t)) + (q* lr (t)qtr(t))) 



lim (q* r (t)qi r (t)) , 

r— *oo 



(H5) 



(q*la(t)qir(t)) 



-27* 



JO 



dtidt 2 R(ti)R*(t 2 ) 



xe 



iz 1 t 1 — lz*t 2 ) 



= e- 2 ^ f dhR 2 c e 2 ^ R2c{1 J 
Jo 2 7 

Substituting this to Eq. (|H4|) . we get 

~ 2 _ 2jk B T 



R n 



(H6) 



(H7) 



APPENDIX Is CALCULATION OF THE 
MOMENTS 

In this appendix we calculate the moments in Eq. 
(Q33J. We have 



d?{q\ - qi *) m {qi - <ti n 0{T)6{T - V) (II) 
= J dY[9\T){ q l q'*r(qi T') 
= J dTKA^LH A\ql qi *) m (qi q[) n }S(r V), 
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where we used the relation ]J H — Lh and L* H = —Lh- 
The quantity inside the brackets in Eq. (fTT|) [which we 
call I] is 

I = (A^LHAHql - q'n m (qi - q[) n 

= -ij^r'e^A^qt - qi *) m (qi " Q'lY |t=o 



ED-ii't-?;) 3 



Z=0 j=0 



(to-0 ! ("-J') !?! j! 



x(- l -)(At)- 1 e ^«*At 9 ™-' g r J lt=o (12) 
Using Eq. (|115|l . we have 



In Eq. I|128|l . by integrating by parts, we can write 

1 k 

dT-^-^G(qi,ql)J2hJkP(T,t) (J2) 



= dr 



dqidql 



G(q 1 ,q* 1 ))Y,bkJke- l(it p^,0)- 



Let us expand 
<9 2 



dqidql 



G(quq* 1 ) = J2 G mnq* 1 m <ll- 



(J3) 



min(m — l.n — I) 

x E 



(m-/)!(n-0! 



a=0 



(to — I — a)!(n — j — a)!a! 



x(e 27t -l) Q |t=o- (13) 

Because of the (e 27 * — 1)° term, the only non-vanishing 
terms in Eq. l"l3|) at i = are for a = or a = 1. So the 
above equation becomes 



).!t).! 



z=o 3=0 



l\j\(m-l)\(n-j)\ 



x(-*|)(e^ t 9l *) m - i (e^ lt (Zl )"^ | t=0 



-1 n-l 



d 



l\j\(m-l-l)\(n-j- 1)! 



x(-i-r )e l{zi - z ^ )t {e l ^ t q* 1 ) m - l - 1 {e- lzlt q 1 ) n - 1 - 1 Y 
x(e 27 '-l)| t=0 

= (-^(e^ql ~ qi ) m (e- tZlt qi - qT |t=o 

+ — Wy(l - e" 27 *) 
at 

x(e^* t a 1 * - a^^-^e-^'o! - a^" 1 | t=0 . (14) 

Substituting Eq. 1(14*1) into Eq. l(TT)) and integrating with 
5(r-r'), we get Eq. (JT23. 



APPENDIX J: FACTORIZATION PROPERTY 

We show the factorization of Eq. (|128f) when p(T, 0) 
has the form 



(Jl) 



We have 



dTq^ql^k-h^PiY^) 

k 

= /dr[( e -^)t g ™ 0l ™E^ J ^*p( r ; o). (J4) 



Since 



L Q J2 h kJk = 0, L v J2 hJk = 0(1/ VZ) (J5) 



and A is expressed in terms of Lq and Ly, 9 — AL#A _1 
treats J2 k bkJk like constant. Neglecting 0(1/ \f~V) terms, 
we can write 



drgr ra a^6 fc J fc e-^p(r,0) 

k 

dr £ bkM(e- l§t )Urq*iTp(r, o). (J6) 



In Eq. """"""J, [(e-^*)tgf gf 1 ]* can be written as (see Eq. 
E"3) 



[(e- m )U? qt n ]* = (^y 1 ( e iLHt A^ql 



*m n 



E 

1=0 



(to - ")!(n- l)W. 



x ^(mz*-nz 1 )t q *m-l q n-l Y l^t ^ ( J? ) 



Since 



(38) 



we can write 



j2hJk[(e- i§t )^Tqrr 

k 

= E WC^^M"]/-* + 0(1/L). (J9) 
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In Eq. I|J9|) . [ ]/_& means that we exclude the fcth field 
mode. With Eq. and neglecting 0(1/ L) terms, Eq. 
(|.T6I) becomes 



dT £ b k J k [(e- M )Urqr}*P^, 0) (J10) 

fc 

= J2 ! dTb k J k [(e- i8t y q ?qr]}- k P(m 

k J 

= J2 I ^6*^[(e-^) t gr«r]/-*5i(Ti)n5*( r *) 
fe ^ fc 

= ^ / dT k b k J k g k (T k ) 

J, J 



fc'^fc 



For any fc we have 



/" dr fc3 (r fe ) = i. 



(Jii) 



Then we can write 

dT k b k J k g k (T k ) 



J dr6 fe j fc5 i(ri) J] 5fc (r fe ) = & fe (J fc ) , (J12) 

dT / _ fe [(e-^)t g ™g 1 *™]/- fc ffl(ri) II ^'( r ^') 

= / dr[(e-^)t<gr]}- fc 9i(ri)n^'( r ^') 

dr^Vi^'p^O), (J13) 



and Eq. <jll0| can be written as 



dTj2hJk[(e- iSt )^TQ*i n rp{r,0) 

k 

Y,h(Jk) J dTql m $e- iSt p{T,Q) 

k 

Y,b k (Jk) [ drqr n q?p(T,t). (J14) 

7.. ^ 



This equation, together with Eqs. I|J2(1 and l|J3(l . leads 
to Eq. (H2EI- 
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